(i-wave pairing from spin fluctuations in the Ka;Fe2_j;Se2 superconductors 
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Angle-resolved photoemission spectroscopy measurements on the recently discovered superconduc- 
tors in the KFe2Se2 family with critical temperatures up to ~ 33K suggest that no Fermi pockets of 
hole character centered on the F point of the Brillouin zone are present, in contrast to all other known 
ferropnictide and ferrochalcogenide superconductors. Using a fluctuation exchange approximation 
and a 5-orbital tight-binding description of the band structure, we calculate the effective pairing 
interaction. We find that the pairing state in this system is most likely to have d-wave symmetry 
due to pair scattering between the remaining electron Fermi pockets at wave vector q ^ {n,n), but 
without any symmetry-imposed nodes for the given Fermi surface. We propose experimental tests 
of this result, including the form of the resonance spectrum probed by inelastic neutron scattering. 

PACS numbers: 74.25.Ha,74.72.-h,75.10. Jm,75.40.Gb 



Introduction. Recently a new family of Fc-chalcogenide 
superconductors Aa;Fe2-ySe2 (A=K,Ca) with Tc oc 30K 
has been discovered These compounds are heavily 
electron doped, such that there are only electron Fermi 
surface pockets, according to angle-resolved photoemis- 
sion (ARPES) studies [2|. The usual argument [3-11 
leading to the most popular "s±" gap structure in the 
Fe-based superconductors requires a F-centered pocket 
to enhance spin fluctuation pairing with wave vector 
Q ~ (tt, 0) in the unfolded (1-Fe) Brillouin zone. In the 
absence of these hole pockets, which were present in the 
previously studied Fe-based superconductors, a gap of 
the s± type is unlikely. Thus the question of the pairing 
mechanism and the structure of the gap in these materials 
remains open. One possibility is that the pairing interac- 
tion is associated with the exchange of spin- fluctuations, 
as considered for the older materials thai the 

effective interaction peaks at a wave vector Q = (tt, tt) 
rather than at (tt, 0). In this case one would expect that 
the gap would have Big (d-wave) symmetry, changing 
sign between the (tt, 0) and (0, tt) electron Fermi surfaces. 
Since there are no portions of the Fermi surface along the 
(tt, tt) direction in the Brillouin zone, there is no symme- 
try reason in a 2D d-wavc gap to have nodes. 

The situation reported by ARPES presents an inter- 
esting new theoretical challenge in these systems, namely 
the calculation of the possible types of superconductivity 
arising in the presence of small pockets of only one type 
of carrier. The proximity of d-wave pairing to the domi- 
nant s± pairing channel in the pnictides generally was 
discussed in earlier spin- fluctuation theories 0, and 
its likelihood in the situation with pockets of one type 
was mentioned briefly in earlier work 0, 0| , but was not 
explored seriously. Recently, Thomale et al.[l3l showed 
that d-wave pairing was likely in the 3K superconductor 



KFC2AS2, which is believed to possess only hole pockets. 
The chalcogcnide analog K^Fe2_j,Se2 is of considerable 
interest not only because the opposite situation obtains, 
but because the critical temperature is an order of mag- 
nitude higher. 



If superconductivity is possible with only one type of 
pocket, there arc in addition two qualitatively different 
situations to address: one in which the hole and elec- 
tron bands overlap in energy, but only one pocket type is 
present due to heavy doping. In addition, a new type of 
situation is suggested by ARPES, one in which an energy 
gap exists between two bands. Doping may then lead to 
a transition between two different symmetry supercon- 
ducting states, or from a superconducting state with one 
symmetry, to an insulating state, and then to another 
symmetry superconductor. These different types of gap 
symmetry transitions as a function of doping will be im- 
portant to explore. 



The main results of our paper are obtained from a fluc- 
tuation exchange RPA calculation for a five orbital tight 
binding model based on an LDA calculation for these Fe 
chalcogen compounds. Within this model, using typical 
interaction parameters, we determine the structure of the 
gap. Wc then discuss ways in which the nodelcss d-wave 
state might be distinguished from ordinary s-wave super- 
conductivity, focussing particularly on the neutron scat- 
tering response that would be expected in such states. 



Model. In the following, we consider a general two- 



FIG. 1. (Color online) The 5-orbital fit to the LDA band 
structure with colors indicating the majority orbital char- 
acter {Ted=dxz, green— dyz, h\ue=dxy, orange=d^2_y2 , and 
magenta=d322_^2 ). The gray points indicate the 10-orbital 
Wannier fit to the full DFT band structure. The splitting 
of the dxz/dyz and the d^cy bands at F has been enlarged to 
remove the hole pockets. 



body onsite interaction Hamiltonian 

i.e i,i'<e 



(1) 



where the interaction parameters C/, C/', J, J' are used 
in the notation of Kuroki et al. 0]. The tight-binding 
Hamiltonian Hq is fitted to the full DFT band structure 
of the parent compound KFe2Se2, calculated within a 
plane wave basis set with ultrasoft pseudopotentials us- 
ing the tools of the quantum espresso package . A 
Wannier projection [12| onto a 10-orbital Fe d basis al- 
lows the determination of the position and the orbital 
composition of the energy bands that can then be fitted 
by a reduced 5-orbital tight-binding model similar to the 
one found for the isostructural BaFe2As2 
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Considering the recent ARPES resuls on K0.8Fe1.7Se2 
reported by T. Qian et al. Q we have artificially en- 
hanced the splitting between the two dxz/dyz bands and 
the two dxy bands at the T point of the backfolded Bril- 
louin zone that is controlled by the nearest neighbor hop- 
ping tx{dxz,dxz) and tx{dxy,dxy), respectively. This al- 
lows us to push the hole pockets below the Fermi level 
without changing the orbital character of the respective 
bands (see Fig. [T]). The effect of a simultaneous down- 
ward shift of the hole bands together with an upward shift 
of the electron bands reducing the difference between the 
bottom of the electron bands and the top of the hole 
bands can be explained if one accounts for the dominant 



FIG. 2. (Color online) The Fermi surface of the 5-orbital 
tight-binding fit with fi — Ef — 0.2. The colors represent the 
major orbital character of the Fermi surface with the same 
color code as in in Fig. [T] 



interband interactions present in the pnictides 1^. We 
have also adjusted the chemical potential to account for 
the reduced electron doping of K2;Fe2-ySe2 with x = 0.8 
and y = 0.3 (0.1 electrons per Fe) compared to the parent 
compound with a: = 1 and y ~ Q (0.5 electrons per Fe). 
In Fig. [2] the Fermi surface for /i = Ep — 0.2 is shown 
with a color encoding of the majority orbital character. 
Note that the square Fermi surface pockets found here 
allow for the possibility of nesting at vectors away from 
(7r,7r). 

To determine the pairing symmetry arising from a 
spin fluctuation exchange picture, we define the following 
scattering vertex r(k, k') in the singlet channel, 



ry(k,k') = Re 



E 



(2) 



X [r,,,,,3,, (k, k', - 0)] ai\ 0^')ai] (-k') 

Here the momenta k and k' are restricted to the electron 
pockets k € and k' G Cj where i and j label either the 
/3i or the P2 Fermi surface, and a'j,(k) arc the orbital-band 
matrix-elements. The orbital vertex functions ^(.^121314, 
represent the particle-particle scattering of electrons in 
orbitals ^1,^4 into i?2,^3 and are given by 



(k,k',c.) 



^C/^xrA(k-k',c.)[/^ 



-[7^xr^(k-k',c.)t/= 



(3) 



The interaction matrices and JJ'^ in orbital space are 
built from linear combinations of the interaction param- 
eters. Their explicit form can be found e.g. in Ref. [§]. 
Here xf^^ and x!^^^ denote the spin-fluctuation con- 
tribution and the orbital-fluctuation contribution to the 
RPA susceptibility, respectively. 

The symmetry function (?(k) of the pairing state can 
then be found by solving an eigenvalue problem of the 
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x'(q) 
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FIG. 3. (Color online) The real part of the RPA enhanced 
static susceptibility for U = 0.96, J — (7/4 in the normal 
state showing a broad peak at q = {n,TT). 



form 



E 



dK 



27rui?(kj|) 



^r.y(k,k')g(k') = Ag(k) (4) 



where the eigenfunction g{\i) corresponding to the largest 
eigenvalue A gives the leading pair instability of the sys- 
tem. 

In the following we parameterize the superconducting 
gap in low-order harmonics and calculate the susceptibil- 
ity in the symmetry-broken state as [isl [l6| 



(5) 



x{G^(fc + q)G''{k) + Ff^i-k - q)F''{k)] 

where the generalized momenta q — (q, aj,„) and k = 
(k, aj„) comprise both the momenta and the Matsubara 
frequencies. The normal and anomalous Green's func- 
tions are given as 

On+U^) A(k) 



G^(fc) 



F^ik) 



(6) 



Here the matrix elements connecting band and orbital 
space determine 

MZuiK q) = «r (k + q)<(k)4'*(k)a;(k + q) (7) 

and the quasiparticlc energies for a band u are given as 
E^(k) — ^yQCk) + A2(k). The inelastic neutron inten- 
sity is proportional to the imaginary part of the spin 
susceptibility in the symmetry-broken state 



(8) 



The multiorbital RPA enhanced spin susceptibility is de- 
fined by a Dyson-like equation 

X^fttiQ,^) = {x'(.q,u;) [l - [7^(9,^)]"'} (9) 

J rstu 
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FIG. 4. The leading d-wave eigenvalue as a function of the 
intra-orbital interaction strength U. Insert: the symmetry 
function g(k) for U = 0.96, J = U/A along the electron pock- 
ets /3i (red) and j32 (blue) as a function of the winding angle 
parameterizing the Fermi vector k relative to the center of 
each pocket. A parameterized fit of g(k) used in the neutron 
scattering calculation is shown as the solid line. 



with interaction matrices If as discussed in Ref. [9|. 

Results. The absence of the hole pocket around F in 
the unfolded 1 Fc per unit cell Brillouin zone removes 
the dominant q = (tt, 0) nesting that is characteristic 
of the LaOFeAs and the BaFc2As2 compounds and is 
responsible not only for the stripe-like SDW instability 
of the undopcd parent compounds but also for the sign 
changing s-wave superconducting state as first pointed 
out by Mazin et al. 0] and later consistently reported 
by RPA, fRG and FLEX calculations for the doped su- 
perconducting materials 0, @, 0, S] ■ With only the elec- 
tron pockets present, the real part of the susceptibility 
reflecting the nesting properties of the electron pockets 
around (tt, 0) and (0,7r) takes a broad plateau-like max- 
imum around q = (tt, tt) that is bordered by two weak 
peaks at q « (tt, 0.6257r) and q w (0.6257r, tt) resulting 
from an enhanced scattering between the flat top and 
bottom parts of the /Si and P2 sheets respectively (com- 
pare Fig.[2|). Here we note that the broad fiat maximum 
corresponding to the nesting of the two electron pockets 
is in contrast to the sharp peak features expected from 
the nesting of an electron and a hypothetical hole pocket. 
The orbital (charge) susceptibility not shown here has no 
pronounced momentum space structure and does not ap- 
proach an instability for the parameters chosen. 

In Fig. |4] we show the leading d-wave eigenvalue as a 
function of U for a Hund's rule coupling J = U /4, a pair 
hopping term J' = J and an inter-orbital Coulomb inter- 
action [/' = [/ — 2 J where the latter two are fixed by the 
local spin-rotational invariance. It reaches an instability 
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around U = 1 corresponding to an anisotropic but node- 
less superconducting gap on the electron sheets as shown 
in Fig. m Due to the lA/mmm symmetry of the crystal, 
the backfolding of the bands from the effective BZ of the 
1 Fe unit cell to the real BZ of the 2 Fe unit cell leads to 
two concentric electron pockets of different size around 
the M point of the backfolded zone. Due to the origin of 
the two concentric electron pockets in the backfolded BZ 
from different positions in the unfolded zone, the super- 
conducting gap exhibits a phase difference of tt between 
them but respects the overall Big symmetry. 

The present calculation was performed using the 1-Fe 
Brillouin zone representation of the band states, which 
neglects the hybridization between the two electron pock- 
ets backfolded onto one another via the body-centered 
cubic symmetry operation translation by (7r,7r,7r). It 
is interesting to ask how the result will change if this 
hybridization or spin-orbit coupling is included. In the 
simplest case, if the unhybridized pockets do not cross 
at fc^ = or ±7r when backfolded, the hybridization will 
affect only the states near ~ ±7r/2, leading to a hor- 
izontal node if the pairing interaction is relati vely in- 
dependent, as we have found in other studies [IJ]. If the 
band structure involves crossings at A:^ = or ±7r, nodes 
with some vertical character may be formed on the elec- 
tron sheets. In contrast to nodes driven by the dominant 
point group symmetry microscopic 2D interaction, how- 
ever, these nodes are a consequence of the crystal space 
group symmetry and their strength is proportional to the 
hybridization between the 1-Fe bands. The density of as- 
sociated quasiparticlc excitations is therefore expected to 
be weak. Further calculations in the 2-Fe zone are un- 
derway to confirm this. 

To calculate the neutron response wc parametrize the 
superconducting gap as A(k) ~ Ao5(k) with 

g{k) = [coskx — cos ky) + 1.62(cos 2/03; — cos2fcy) (10) 

This fit is shown as the solid line in the inset of Fig. SI In 
Fig. [5] we show the imaginary part of the RPA enhanced 
susceptibility in the d-wave state for a momentum trans- 
fer of q = (7r,0.6257r) in (a) and q = (tTjTt) in (b) in 
the IFe/unit cell Brillouin zone. Here the low energy 
spectral weight is suppressed upon the opening of the 
superconducting gap and is transferred to higher ener- 
gies. The pronounced resonance peak around w = 2A 
appears only for a relative sign change of the gap on 
the two electron pockets such that the coherence factor 
1 - A(k)A(k + q)/(£'(k)£;(k + q)) does not vanish. 

There are proposals that an s-wave gap may arise from 
the orbital term in Eq. [3] when local Fe phonon modes 
are included [13, In Fig. [5] we have added results 

for an s-wave gap taken equal to the average magnitude 
of the d-wave gap for comparison. While an s-wave gap 
could well have anisotropic structure, we expect that the 
difference in between a Big and an Aig gap will 

be significant. In particular if the orbital fluctuations 
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FIG. 5. (Color online) The imaginary part of the RPA-BCS 
dynamic spin susceptibility versus lj for q — (7r,0.6257r) (a) 
and q = (tt, tt) (b) for the normal, d-wave and s-wave states. 
The interaction parameters were chosen as U = 0.96 and J — 
U/4. 

are dominant, the response in the magnetic scattering 
channel will be further suppressed compared to the d- 
wavc response illustrated in Fig. [S] 

Conclusions. We have argued, based on an RPA 
treatment of a generalized multiorbital Hubbard model, 
that the absence of the F-ccntercd hole pocket in the 
KFe2Se2 superconducting materials should lead directly 
to a strong d-wave pairing instability without nodes on 
the remaining M-centered electron pockets. The appear- 
ance of d-wave pairing in this family of unconventional 
superconductors in the limit when only one pocket is 
present would be strong support for pairing by spin fluc- 
tuations in these systems. It appears to us that the mea- 
surement of a peak in the inelastic neutron scattering 
spectrum near tt, tt would be the easiest way to test this 
prediction. 

Since the inelastic neutron scattering is mostly sensi- 
tive to the Fe lattice, it is possible to distinguish with this 
technique between the low q and the q = (tt, tt) scatter- 
ing in the unfolded IFc/unit cell Brillouin zone, although 
both signals would be backfolded on the F point in the 
2Fc/unit cell relevant e.g. for the interpretation of the 
angle resolved photoemission results. Therefore the pro- 
posed experiment can provide a direct measurement of 
the q vector dominating the repulsive interaction and 
eventually leading to a sign change of the superconduct- 
ing gap on Fermi surface regions connected by q. 

In a recent preprint using a functional renormalization 
group approach, F. Wang et al. also conclude that the 
leading pairing instability of a Ka:Fe2-j/Se2 model occurs 
in the dx2-y2 channel [l9[. 
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